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Abstract 

New calculations to over ten million time steps have revealed a 
more complex diffusive behavior than previously reported, of a point 
particle on a square and triangular lattice randomly occupied by mir- 
ror or rotator scatterers. For the square lattice fully occupied by 
mirrors where extended closed particle orbits occur, anomalous dif- 
fusion was still found. However, for a not fully occupied lattice the 
super diffusion, first noticed by Owczarek and Prellberg for a particu- 
lar concentration, obtains for all concentrations. For the square lattice 
occupied by rotators and the triangular lattice occupied by mirrors or 
rotators, an absence of diffusion (trapping) was found for all concen- 
trations, except on critical lines, where anomalous diffusion (extended 
closed orbits) occurs and hyperscaling holds for all closed orbits with 

universal exponents df = — and r = — . Only one point on these 

critical lines can be related to a corresponding percolation problem. 
The questions arise therefore whether the other critical points can be 
mapped onto a new percolation-like problem, and of the dynamical 
significance of hyperscaling. 



KEY WORDS: Diffusion, Lorentz lattice gas, critical point, hyper- 
scaling, super diffusion. 
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1 Introduction 

In a number of previous papers the diffusive behavior of Lorentz Lattice Gas 
Cellular Automata (LLGCA) has been studiedt^"^^. Here a point particle moves 
in (constant) discrete time steps on a discrete lattice from site to site, a number 
of which is occupied randomly by stationary scatterers which scatter the particle 
according to strictly deterministic rules. The nature of the diffusive process of the 
particle through the scatterers has been the object of these investigations. In this 
paper we will confine ourselves to the motion of a particle on a square or triangular 
lattice occupied by fixed scatterers which remain unchanged during the diffusion. We 
will consider two different scatterer models: the mirror model and the rotator model. 
In the mirror (rotator) model a particle is scattered upon collision with a mirror 
(rotator) to the right or the left, depending on whether the mirror (rotator) is a right 
or a left mirror (rotator), respectively (cf.fig.I). The fraction of right (left) scatterers 
(mirrors or rotators, respectively) on the lattice will be denoted by Cr (Cl), so that 
Cr + Cl = C is the total fraction of lattice sites occupied by scatterers, i.e., the 
concentration C of the scatterers. We also choose the time step = the lattice distance 
= speed = I. For a given random placement of the mirrors (rotators), a particle will 
describe - from a given initial position - a random-like walk through the lattice (see 
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fig.l). The diffusive behavior of the particie wiii be obtained by averaging over aii 
possibie random configurations of the mirrors (rotators). It can be characterized by 
a number of quantities, of which we will only consider the mean square displacement 
A(t) and the radial distribution function P(r, t), defined by: 



respectively. Here r = (r, 6) is the position of the particle in polar coordinates with 
respect to the origin, so that r = |r| is the distance of the particle from the origin. 
The average in eq. (1) is over all random configurations of the scatterers at fixed 
Cr and Cl and the sum in eq.(2) over all possible angles consistent with the lattice. 
P[f^t) is the probability to find the particle at the position r at time t, given that it 
was at the origin at time t = 0, so that P(r, t) is the probability to find the particle 
a distance r away from the origin at time t. 

From eq. (1), we can define a time dependent diffusion coefficient Dit), by: 



so that a diffusion coefficient D exists if in 



A(t) =< r\t) > 






D = lim D{t) 



(4) 



3 



the right hand side (r.h.s.) exists. 

When the diffusion is normai (we wiii caff it ciass I), P(r, t) wiii be a Gaussian 
given by: 

f r2 

PMrA) = — ^=exp( ) (5) 

which impiies A(t) ~ t, but not vice versa (cf. tabic I). 



2 The Mirror Model 

For the mirror modei - which was introduced by Ruijgrok; and one of us^^^ - 
it was found that, although D existed, P[f^t) was not Gaussian. This was called 
anomalous (or class II) diffusion^ (cf.figs.2a-b). This non-Gaussian behavior is due 
to the presence of closed (periodic) orbits, especially near the origin. This is in 
contrast to normal (Gaussian) diffusion where there are no closed orbits, but only 
open trajectories of the particles. The time evolution of P is indicated in fig. 2b; it is 
dominated by closed orbits and zig-zag motions (cf.fig.2c). Since the time at which 
orbits close can be arbitrarily large, (infinite) extended closed orbits can occur. 

The question arises whether all trajectories eventually close for all C and Cr/Cl- 
For C = 1 this question was answered affirmatively by a theorem of Bunimovich and 
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Troubetzkoyt'']: if C = 1 and Cr > 0; C'l > then all trajectories are closed with 
probability 1. In this case (C = 1) one even knows how fast the trajectories close 
since it has been shown that the probability Po{t) to find an open orbit after time t 
IS given 

Po{t) ~ r'- (6) 

so that the number of open orbits decreases in time according to eq.(6) and all trajec- 
tories eventually close. The result (6) was obtained by noting a connection between 
the dynamical problem considered here and a bond percolation problem on the square 
lattice (see below). 

For C < 1, only numerical evidence existst^^ that 

^"") Hik 

where Ci and C2 are constants. (7) implies that all trajectories still close, albeit much 
slower than (6), as illustrated in fig. 3 (cf.ref. [9]). 

Earlier work suggested anomalous diffusion for all C, which leads to a dynamical 
phase diagram of the mirror model on the square lattice like that in fig. 4a. This phase 
diagram was based on calculations extending typically to 4,000 - 10,000 time steps. 
The method of calculation is described in refs. [2-5] and in the few calculations done 
for a larger number of time steps the class II behavior was consistent with the error 
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bars of the data points (cf.fig.2a)^. 

However, in a recent paper by Owczarek and Prellberg[^°] it was shown that for 

2 1 . . 

the particular case ^ ~ ^ ~ = ~, the diffusion process was not anomalous 

but super diffusive, since Dit) grew logarithmically with t. This result was based on 
very (many months) long calculations of up to a million time steps, giving results 
with very small error bars. This prompted us to extend all our calculations to longer 
times, in order to see how prevalent this super diffusive behavior was and to correct 
our phase diagram. In order to reduce the calculation times and in addition avoid 
the use of periodic boundary conditions, we employed a modification of the method 
developed by Ziff, Cummings and Stellt^^l. Here a virtual lattice of 65536 X 65536 
sites is divided into small blocks (256 X 256 lattice sites each) where scatterers are 
initially put randomly only on that block where the particle starts its motion, while 
no scatterers will be put on all the other blocks until the particle enters. This provides 
an enormous saving in memory and (along with other efficiencies) allows a reduction 
in computer time from months to days for obtaining results with very small error bars 
without any boundary effects. The details will be given in the next paper^^^l. 

Figs.5a-b show our results for Dit) for C = f as well as the decrease in the number 

special case arises when Cr (or Cl)=0 and < Cl (or Cr)< 1, and (zig-zag) 
propagation occurs along the direction of the mirrors, while Gaussian (in fact, Boltzmann) 
diffusion takes place along the direction perpendicular to the mirrors. 
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No{t) of open orbits for 10,000 particles for a number of values of Cr/Cl up to a 
million time steps. The statistical errors were determined by doing the calculations 
in two steps: first an average was made over all 10,000 particles, with a different 
random configuration of the scatterers for each particle, then further averages were 
computed over typically three runs, involving three samples of 10,000 particles for 
each. The standard deviations of the mean are plotted as the error bars of the data 
in the figures. If the error bar does not appear, the error bar is inside the symbol. 
These results are consistent with class II behavior. Fig. 6a shows D[t) for a number 
of concentrations for C < 1 and Cr = Cl- contrary to what was found in fig. 2a, D[t) 
increases for sufficiently long times logarithmically in time for all C < 1, according 
to: 



where A is essentially independent of Cr for all < C < 1 (cf.fig.6b) (for C very close 
to or 1, we need to run for much longer times than feasible to obtain a constant 
slope). Fig. 6c shows D[t) for a number of concentrations Cr/Cl for C = 0.8. The 
logarithmic increase (8) of D[t) is clearly visible, but for Cr ^ Cl, the constant A 
is Ci?-dependent (cf.fig.6d). Thus, while the coefficient A appears to be universal 
for Cr = Cl it is not for Cr ^ Cl- As a result of the behavior of Dit), the phase 



D{t) ~ Alogt 
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diagram for mirrors on a square lattice is not that of fig. 4a, but that shown in fig. 4b. 

The P(r, t) for C = 1 and Cl = Cr = 0.5 is plotted in fig. 7a, while that for 
C = 0.6 and Cl = Cr = 0.3 is plotted in fig. 7b for a number of different times^. In 
both cases the curves increase in length and decrease in height with increasing time. 
Although the first case is that of anomalous diffusion (class II) and the second case of 
super-diffusion, the difference between the two curves is only significant around r 0, 
where in the first case almost one order of magnitude more closed orbits occur than 
in the second case. For larger times, a difference is hardly noticeable and it is not 
clear whether this is due to the difference in concentration or in diffusive behavior. 

We believe that the origin of the super diffusive behavior for C < I is two-fold: 

1. the slow decay of the number of open trajectories, as given by eq.(7); 

2. the possibility of large zig-zag motions (cf.fig.2c), because of the presence of unoc- 
cupied sites on the lattice. 

3 Mirror Model and Percolation 

For later discussions, we now sketch how the trajectories of the particle on the 

fully occupied square lattice (i.e., for C = I) in the mirror model can be related to 

^There is no obvious difference in behavior from fig. 2b, except that we used periodic 
boundary conditions in fig. 2b, rather than the procedure of ref.[ll]. 
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percolation clusters of a bond percolation problem on two sublattices of the original 
lattice[2'i3]. In fig.8 the square lattice with two square sublattices is shown. We note 
that each lattice site of the original lattice is part of both sublattices"^. Since only one 
mirror can be put at a lattice site, the mirror at this site can belong to only one of the 
two ambient sublattices. As a result, each sublattice is only half filled with mirrors, 
so that the probability p that a bond in each sublattice is occupied by a mirror is 
— . This is just the critical value for bond percolation on a square lattice, so that the 
mirrors are at the bond percolation threshold on the two sublattices. Thus, if the size 
of the mirrors is chosen equal to the bond length of the sublattices, bond percolation 
clusters appear. We note that each closed orbit with period larger than 4 time steps 
is formed by refiection between the inner and outer perimeters of bond clusters on the 
sublattices (cf.fig.8). There is no ambiguity about "inner" or "outer", as the closed 
orbits do not cross themselves and the inner and outer perimeters refer to clusters on 
different sublattices. This mapping implies the following analogy. From percolation 
theory we know that the mean square gyration radius < i?^ >o of the perimeter of a 
percolation cluster with N perimeter sites grows with N as: 

< R% >„~ = m (9) 

^The sublattices were plotted wrongly in fig.l in ref.[2]. 
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where the fractal dimension df = —. Furthermore the probability Po{N) to find an 
open percolation cluster with N perimeter sites decreases with N as: 

Po{N) ~ iV^-^ = N-'- (10) 

15 

so that the size distribution parameter r = — . We note that dj and r satisfy a 
hyperscaling relation: 

Similarly, for the trajectories of the particles on the lattice, one finds that the mean 
square displacement Ao{t) for particles on open trajectories at time t increases as: 

A„(t)~tt (12) 

while the probability for an open trajectory at time t is given by: 

Po{t) ~ t-'- (6) 

The analogy of (9), (12) and also (10), (6), respectively, is obvious and illustrated in 
fig. 5b. 
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4 The Rotator Model 

We now discuss our rotator model, which is a special case ol a set ol rotator models 
introduced by Gunn and Ortuno^^'^^. The scattering rules are illustrated in fig.l and 
although similar to those of the mirror model, differ from those in their scattering 
results in half of the cases (cf. the rules in fig.l); they do not lead to time-reversible 
particle motions, as in the mirror model. 

In earlier work, the rotator model seemed to behave for Cr = Cl very similarly 
to the mirror model (cf.fig.9), while for Cr ^ Cl a qualitative difference was no- 
ticed. Then, instead of showing a class II anomalous diffusive behavior as the mirror 
model did, a dynamical phase transition was observed, for sufficiently large Cr or 
Cl, from anomalous to an absence of diffusion-to which we will refer as no-diffusion- 
where the mean square displacement became bounded (A(t) < constant) and all 
particles trapped, so that no extended closed orbits occurred. We called this class IV 
behavior^. This behavior is illustrated in fig. 10a for D[t) and in fig. 10b for P[r^t). 
We note, that indeed for t > 2^°, the curve for D[t) on logio-logio scale has a slope 
of -1 and P appears stationary and does not seem to change anymore, as all parti- 
cles are trapped by that time. This leads to a phase diagram as pictured in fig. 11a, 
where the phase transition occurs across two (approximately determined) straight 
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lines, anchored on the Cr and Cl axes at Cr = Cl = 0.593, the critical concentra- 
tion for site percolation on the square lattice. The point Cr = Cl = 0.5 for C = 1, 
which can be mapped on the corresponding point for the mirror model, has class II 
behaviort^l. This result was consistent with another theorem proved by Bunimovich 
and Troubetzkoyt''] which stated that 1) there exists a critical concentration Cr^^ or 
Clc^ ^ (0-5, 1) (our 0.593), such that, for Cr > Cr^^ or Cl > Clct^ trajectories are 
periodic with probability 1; 2) for C = 1 all trajectories are periodic with probability 
1. The proof is based on the observation that if an infinite percolation cluster of, 
say, right rotators exists, the particle inside the cluster is everywhere surrounded by 
a contour consisting of only one kind of rotators and consequently is trapped. The 
regions covered by this theorem are indicated in fig. 11a. 

Recently we have extended all calculations for this model also to a million and 
sometimes close to more than ten million time steps and again found very different 
results as were obtained before for typically 4,000 - 10,000 time steps. 

In fig. 12 Dit) is plotted for C = 1.0 and 0.8 for Cr = Cl as well as for a particular 
value Cr ^ Cl for C = 0.9. While for Cl = Cr and C = 0.8, D{t) 0, i.e., class 
IV behavior or trapping is observed, for the particular concentrations Cr = 0.423, 
Cl = 0.477 for C = 0.9, D[t) seems to approach a constant ^ 0, i.e., class II behavior 
or extended closed orbits occur, like for Cl = Cr = 0.5 when C = 1. The possible 
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occurrence of extended closed orbits for special values of Cr and Cl is confirmed 
by further study. Fig. 13a shows the number of open orbits No{t) for C = 1 for 
a variety of Cr and Cl- Only for Cr = Cl = 0.5 do extended trajectories occur 
and a corresponding slow decay of No{t) occurs, while for Cr ^ Cl, the precipitous 
decay of No{t) clearly indicates the sudden trapping of particles when t approaches 
a critical value. Similarly figs.l3b-c show class II behavior only for the particular 
values Cr = 0.455, Cl = 0.495 and Cr = 0.46, Cl = 0.39 for C = 0.95 and C = 0.85, 
respectively. The number of closed orbits Nc{t) at these critical concentrations of 
Cr and Cl is a minimum when compared to those at all other concentrations of Cr 
and Cl at a given value of C. For all these cases, the decay of No{t) with time is 
again given by eq.(6), confirming the suspected class II behavior deduced from the 
behavior of D[t). Fig. 13d shows the similarity of the decay of No{t) for a number of 
concentrations and suggests a dynamical phase diagram as given in fig. lib, instead 
of that in fig. 11a. The existence of critical points for C < 1 was noticed before by 
Ortuno, Ruiz and Gunnt^^^, but the location and properties of their points are quite 
different from ours. 

The present situation can be summarized as follows. 
1. There appear to be two symmetric critical lines - when Cr/Cl is critical, so is 
Cl/Cr - in the phase diagram, where extended closed orbits of a particle can occur. 
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Outside these critical lines there is trapping everywhere. 

2. Only for C = 1 is there a connection of the rotator model with percolation through 
a mapping on a bond percolation problem^^'*]. 

In fact, at C = 1 a given closed orbit can always be mapped into the perimeter 
of a bond cluster as follows: we replace each rotator along the trajectory by a bond 
which looks like a mirror with a length of \/2 lattice distances and with an orientation 
such that the bond does not cross the particle trajectory (cf.fig.8). Like in the mirror 
model, the particle trajectory bounces back and forth between the inner and outer 
perimeters of bond clusters. Like for the mirror model, there is no ambiguity about 
"inner" and "outer" perimeters since the closed orbit does not cross itself and the 
bonds do not cross the closed orbits either and the inner and outer perimeters of 
bond clusters reside on different sublattices, respectively (cf.fig.8). However, unlike 
for the mirrors, the inner perimeter of a bond cluster always corresponds to one kind 
of rotators (R or L) while the outer perimeter of bond cluster always corresponds to 
the other kind of rotators (L or R). 

We will now use this to argue that for C = 1 percolation occurs only for Cl = 
Cr = 0.5, unlike for the mirrors. First we note that the number of (the same kind of) 
rotators corresponding to the perimeter of the inner bond cluster is always smaller 
than that of the corresponding outer bond cluster. Because the particle trajectory 
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has a one to one relation to the inner bond cluster, it follows that for large closed 
orbits the probability for a bond of the inner perimeter of bond cluster will depend 
on the smaller of the two concentrations Cr and Cl at C = 1. Therefore, since 
Cl = Cr = 0.5 just corresponds to the bond percolation cluster threshold, where the 
probability p = Pc = 0.5, so that there will exist extended closed orbits, leading to a 
diffusive behavior of class II. However, for Cl 7^ Cr^ the perimeter of the inner bond 
cluster corresponds to a p < pc = 0.5, so that there is then no inner bond percolation 
cluster and there are no extended closed orbits, i.e., all particles are trapped and the 
diffusive behavior is that of class IV. 

For C < 1, no such mapping seems to exist, since no percolation clusters can 
be formed because of the empty sites on the lattice. Nevertheless, we determined 
independently that r = y by Po{t) ~ t'^~^ = t~~ (cf.figs.l3a-d) as well as = | 
from Po{t)Ao{t)/t ~ t'rf/f /t = t~rtr /t ~ constant (cf.fig.l4a) leading, within the 
experimental errors, to the validity of the hyperscaling relation eq.(ll) along the two 
symmetric critical lines. 

3. In spite of what one would expect physically, for C < 1, there appears to be no 
obvious influence of the site percolation transition and the formation of an infinite 
cluster of right (or left) rotators, on which the proof of the above mentioned Buni- 
movich - Troubetzkoy theorem is based. The Bunimovich - Troubetzkoy theorem 
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seems therefore only sufficient, not necessary for trapping of the particle (class IV 
behavior) to occur. 



5 The Mirror and Rotator Model on the Trian- 
gular Lattice 

A piece of a triangular lattice is shown in fig. 15a. The motion of the particle takes 
place on a triangular lattice if the scattering of the particle occurs over the largest 
possible angle 27r/3, but on a honeycomb (sub) lattice, if the scattering occurs over 
the smallest angle 7r/3t"^^. The latter case will be discussed in the following paper^^^l. 
Right and left mirrors as well as rotators can be defined on the triangular lattice as 
illustrated in fig. 15b. 

For relatively short times [t < 10,000), we found for both (mirror and rotator) 
models the dynamical phase diagram of fig. 16a, similar to that for rotators on the 
square lattice (cf.fig.lla). The only difference is that the critical concentration for site 
percolation on the triangular lattice is Cr^^ = Clct = 0-5 instead of 0.593. A similar 
theorem as for the rotator model on the square lattice has been proved for the mirror 
and rotator models on the triangular lattice by Bunimovich and Troubetzkoyt'']. The 
regions covered by this theorem are indicated in fig. 16a. 

The recently carried out extended calculations up to a million or more time steps 
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again showed very different resuits. 

In fig. 17a, D[t) is plotted for a number of values of Cl and Cr sliowing tliat, 
for Cl = Cr^ Dit) approaclies a constant, wliile for Cl 7^ Cr it decays witli a 
slope of —1, indicating, with eq.(3), that the mean square displacement is bounded. 
There is again the occurrence of a critical concentration at Cl = Cr for C = 0.8 as 
illustrated in fig. 17b, where again the precipitous decay of No{t) for Cl 7^ Cr clearly 
indicates the sudden trapping of particles for a finite t. The slow power-law decay of 
No{t) according to eq.(6) for Cr = Cl = 0.4 is due to the occurrence of (infinitely) 
extended closed orbits and holds for all Cl = Cr. Fig. 17c shows this behavior for a 
number of different C for Cr = Cl and fig. 16b gives the new phase diagram, which 
is - for the triangular lattice - identical for the mirror and the rotator models. 

This same behavior of mirror and rotator model can be understood by replacing 
all the right (left) rotators on a particle trajectory by either right (left) mirrors or left 
(right) mirrors, while not changing the empty sites on the particle trajectory. Thus 
in fig. 18, we only need to replace all right mirrors by right (left) rotators and all left 
mirrors by left (right) rotators, to obtain the same trajectory with the same (opposite) 
direction of particle motion, as in fig. 18, respectively. So, the mapping will be either 

^rotator ^mirror ^rotator ^mirror ^rotator ^mirror ^rotator ^mirror ^hcrC 

L L ^ R R L R ^ R L ' 

is no difference between these two cases since the diffusive behavior of the models is 
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invariant for an interchange of Cl and Cr. 

We notice in fig. 16b that there is now only one critical line Cr = Cl, which 
exhibits class II behavior. Like for the square lattice we determined independently 
the exponents r = y and df = ^ (cf.figs.l4b and 17b-c, respectively) leading, within 
the experimental errors, to a verification of the hyperscaling relation eq.(ll) along the 
entire critical line. For all other concentrations class IV behavior obtains. Only for 
C = 1 for both (mirror and rotator) models is there a mapping to a site percolation 
problem on the same lattice obtained by connecting the same kind of nearest neighbor 
scatterers (mirrors or rotators) (cf.fig.18). 

The present situation can be summarized as follows. 

1. There appears to be one critical line in the phase diagram where (infinite) extended 
closed orbits of a particle occur and hyperscaling relation eq.(ll) holds. Outside this 
critical line there is everywhere trapping. 

2. Only for C = 1 for both (mirror and rotator) models is there a connection with 
site percolation on the same lattice as that of the particle trajectory (fig. 18) (cf. in 
section 4, the mapping of the particle trajectories for the rotator model on the square 
lattice onto bond clusters on dual lattices). 

3. For C < 1, there appears no noticeable infiuence of the percolation transition, 
on which the above mentioned proof of Bunimovich-Troubetzkoy theorem is based. 
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Thus, the Bunimovich-Troubetzkoy theorem seems also in this case only a sufficient, 
not a necessary condition, for class IV behavior to occur. 

6 Outlook 

We summarize our new results in table II. 

We close with the following questions and remarks. 

1. The most striking result of these investigations seems to be the existence for both 
lattices of critical lines which are extentions of percolation transitions and appear 
to have universal percolation perimeter critical exponents satisfying the same hyper- 
scaling relation eq.(ll) as found for the percolation problem. On these critical lines 
extended closed orbits occur, but what determines for a given concentration C the 
critical concentrations for these extended orbits to occur, is unknown to us. 

2. Are there for a not fully occupied lattice (C < 1) generalized percolation cluster- 
like structures to which the extended closed orbits are geometrically related, like the 
percolation clusters on a fully occupied lattice (C = 1), i.e., does our dynamical model 
suggest something new - a possible generalization - for the percolation problem? 

3. What is the dynamical meaning of hyperscaling relations for non-percolation re- 
lated trajectories? 
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4. Can one understand that in all cases -, i.e., on both the square and the triangular 
lattice and for both mirror and rotator models (except for the mirror model on the 
square lattice for Cr or Cl = 0) - the trajectory of a moving particle eventually 
always closes? Furthermore, how can one see that in most cases the scatterers will 
lead to a quick trapping of the particle, but that in some - rather special - cases 
closed orbits of any extent can occur so that the closing of orbits proceeds power-law 
or logarithmically slow. We note that the Bunimovich-Troubetzkoy theorems proved 
so far do not distinguish between class II and class IV closing of orbits. 

5. In all cases considered here, all trajectories close eventually. This does not imply 
that no-diffusion takes place in any form; for anomalous (class II) diffusion, one can 
still define a finite diffusion coefficient by the relation (4). In fact, a whole range of 
possible diffusive behaviors occurs ranging from super-diffusion to no-diffusion. The 
closing of all trajectories does therefore not by itself say anything about the way 
the particles diffuse - as measured by their mean square displacement - through the 
scatterers. 

6. The previously obtained results for the hipping mirror and rotator models^^l do 
not seem to be affected by the present extension of the time scale of our calculations. 
They remain therefore unchanged. 
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Figure Captions 

Fig.l Typical mirror (a) and rotator (b) configurations, particle trajectories and scat- 
tering rules for the mirror and rotator models, respectively, on the square lattice. 
Fig. 2 (a) Old calculation of the diffusion coefficient D as a function of the time t on a 
logio-logio scale for the fixed mirror model on the square lattice for Cl = Cr = 0.25^^^; 
(b) corresponding radial density distribution functions P(r, t) as a function of distance 
r from the origin at time steps t = 2^° (O), t = 2^^ (+) and t = 2^"^ (□), respectively. 
The stationary peak near the origin is due to closed orbits; (c) a possible closed orbit 
(A) and zig-zag motion (B) for the mirror model on the square lattice for C < 1. 
Fig. 3 Inverse probability for open orbits as a function of log2t for the mirror model 
on the square lattice for Cl = Cr = 0.4. 

Fig. 4 (a) Old phase diagram for the mirror model on the square latticet^^; (b) new 
phase diagram. 

Fig. 5 (a) Diffusion coefficient D as a function of logiot for the mirror model on the 
square lattice for C = 1 at Cl = Cr = 0.5 (O), Cl = 0.55, Cr = 0.45 {+) and 
Cl = 0.75, Cr = 0.25 (□) (the first two overlap); (b) number of closed orbits as a 
function of t on a log2-log2 scale for the mirror model on the square lattice for C = 1 
at Cl = 0.6, Cr = 0.4 (O) and Cl = 0.5, Cr = 0.5 (+). The fines through the points 
are drawn to guide the eye, virtually coincide. The slope for both curves —1/7. 
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Fig. 6 (a) Diffusion coefficient D as a function of iogiot for the mirror modei on the 
square iattice for Cl = Cr = 0.4 (O), Cl = Cr = 0.35 (+) and Cl = Cr = 0.3 (□); 

(b) siope of D/logiot as a function of the concentration of mirrors C for Cl = Cr] 

(c) as in (a) for Cl = Cr = 0.4 (O), Cl = 0.5, Cr = 0.3 (+), Cl = 0.6, Cr = 0.2 
(□), Cl = 0.7, Cr = 0.1 (x) and Cl = 0.75, Cr = 0.05 (A); (d) slope of D/logwt 
as a function of concentration of left mirrors Cl for C = 0.8 and Cl 7^ Cr. The lines 
through the points are drawn to guide the eye. 

Fig. 7 Radial distribution functions P(r, t) for the mirror model on the square lattice 
as a function of the distance r from the origin at time steps t = 2^^ (O), t = 2^"* (+) 
and t = 2^^ (□), respectively, (a) for C = 1 at Cl = Cr = 0.5 (peak value of about 
0.27 is not shown in the figure); (b) for C = 0.6 at Cl = Cr = 0.3. 
Fig 8. A typical particle trajectory (thin solid line with arrows) on the square lat- 
tice (thin solid lines) relates to the perimeters of bond clusters each on one of two 
sublattices, respectively (dashed and dotted lines, respectively). The thick lines are 
the perimeters of the bond clusters (also the mirrors for the mirror model), while "L" 
and "R" refer to rotators for the corresponding situation for the rotator model. The 
particle trajectory resides between the inner and outer perimeters of the two clusters. 
Fig. 9 Old calculation of diffusion coefficient D as a function of t on a logio-logio scale 
for the rotator model on the square lattice for Cl = Cr = 0.25^. 
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Fig. 10 (a) Diffusion coefficient D as a function of the time t on a iogio-iogio scaie for 
the rotator modei on the square iattice for Cl = 0.6, Cr = 0.2; (b) corresponding 
radiai distribution functions P(r, t) as a function of distance r from the origin at time 
steps t = 2^ (O), t = 2^° (+) and t = 2^^ (n)^ respectively. 

Fig. 11 (a) Old phase diagram for the rotator model on the square lattice^; the part 
with dashed lines corresponds to that covered by the theorem of Bunimovich and 
Troubetzkoy; (b) new phase diagram. The dash-dotted lines indicate the extrapo- 
lated behavior we expect for C < 0.7, the minimum value of C we considered. 
Fig. 12 Diffusion coefficient D as a function of t on a logio-logio scale for the rotator 
model on the square lattice for the critical concentrations for C = 0.9 at Cl = 0.477, 
Cr = 0.423 (O), C = 1 at Cl = Cr = 0.5 (+) and C = 0.8 at Cl = Cr = 0.4 (□); 
critical concentrations are printed in bold face. 

Fig. 13 Number of open orbits out of 10,000 trajectories as a function of t on a log2- 
log2 scale for the rotator model on the square lattice, (a) C = 1 at Cl = Cr = 0.5 
(O), Cl = 0.51, Cr = 0.49 (+), Cl = 0.52, Cr = 0.48 (□), Cl = 0.53, Cr = 0.47 
(x) and Cl = 0.54, Cr = 0.46 (A) (the number of closed orbits has a minimum at 
Cl = Cr = 0.5); (b) C = 0.95 at Cl = Cr = 0.475 (O), Cl = 0.485, Cr = 0.465 
(+), Cl = 0.495, Cr = 0.455 (□) and Cl = 0.505, Cr = 0.445 (x) (the number 
of closed orbits has a minimum at Cl = 0.495, Cr = 0.455); (c) Cl = Cr = 0.425 
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(O), Cl = 0.43, Cr = 0.42 (+), Cl = 0.44, Cr = 0.41 (□), Cl = 0.45, Cr = 0.4 
(x), Cl = 0.46, Cr = 0.39 (A) and Cl = 0.47, Cr = 0.38 (*) (the number of 
closed orbits has a minimum at Cl = 0.46, Cj? = 0.39); (d) same as in (c) for 
Cl = Cr = 0.5 (O), Cl = 0.495, Cr = 0.455 (+), Cl = 0.477, Cr = 0.423 (A), 
Cl = 0.46, Cr = 0.39 (□) and Cl = 0.44, Cr = 0.36 (x). The lines through the 
points are drawn to guide the eye. 

Fig. 14 Contribution of diffusion coefficient from open orbits Po{t)Ao{t)/t as a func- 
tion of time t on a log2-log2 scale for (a) C = 0.85 at Cl = 0.46, Cr = 0.39 (O), 
Cl = 0.45, Cr = 0.4 (+) and Cl = 0.47, Cr = 0.38 (□) for the rotator model on the 
square lattice; (b) C = 0.8 at Cl = Cr = 0.4 (O) and Cl = 0.41, Cr = 0.39 (+) 
for the rotator model on the triangular lattice. The approach to a horizontal line is 
apparent for the critical concentrations. 

Fig. 15 (a) A typical particle trajectory on the triangular lattice, three-sided dotted 
lines stand for right and left mirrors and "R" and "L" stand for right and left rota- 
tors, respectively; (b) the scattering rules for the mirror and rotator models on the 
triangular lattice. 

Fig. 16 (a) Old phase diagram for both mirror and rotator model on the triangular 
lattice, dashed lines as in fig. 12(a); (b) new phase diagram. 

Fig. 17 (a) Diffusion coefficient D as a function of t on a logio-logio scale for the ro- 
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tator model on the triangular lattice for C = 1 at Cl = Cr = 0.5 (A), Cl = 0.45, 
Cr = 0.55 (x), C = 0.8 at Cl = Cr = 0.4 (*), C = 0.79 at C'l = 0.4, Cr = 0.39 
(O), C = 0.65 at Cl = Cr = 0.325 (+) and C = 0.5 at Cl = Cr = 0.25 (□); (b) 
number of open orbits out of 10,000 trajectories as a function of t on a log2-log2 
scale for the rotator model on the triangular lattice for C = 0.8 at Cl = Cr = 0.4 
(O), Cl = 0.41, Cr = 0.39 (+), Cl = 0.43, Cr = 0.37 (□), Cl = 0.44, Cr = 0.36 
(x) and Cl = 0.45, Cr = 0.35 (A) (the number of closed orbits has a minimum at 
Cl = Cr = 0.4); (c) same as in (b) for Cl = Cr = 0.5 (O), Cl = Cr = 0.4 (+), 
Cl = Cr = 0.3 (□) and Cl = Cr = 0.25 (x). The lines through the points are 
drawn to guide the eye. 

Fig. 18 A typical particle trajectory (closed orbit) on the triangular lattice (thin and 
thick solid lines with arrows) relates to site clusters on the same triangular lattice 
as for the particle trajectory. The thick lines are the perimeters of the site clusters. 
The three-sided dotted lines represent the mirrors for the mirror model, "L" and "R" 
represent the rotators for the rotator model and the particle trajectory is between 
the inner and outer perimeters of the clusters. We can see that the mirror model and 
the rotator model are the same by replacing the right (left) mirrors by right (left) 
rotators or vice versa. 
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Table Captions 

Table I - Comparison of normal and anomalous diffusion 
Table II - Comparison of different diffusive behavior 
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Table I — Comparison of normal and anomalous diffusion 





Normal 
(class I) 


Anomalous 
(class II) 


Ait) 


~ t 


~ t 




Gaussian 


Non-Gaussian 
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Table II — Comparison of different diffusive behavior 



Lattice 


MIRROR 


ROTATOR 


Square 


< C < 1 — > Super-diffusion 


< C < 1 ^ No-diffusion (class IV) 
except for 2 critical lines — > 
Anomalous diffusion (class II) 


C = 1 — > Anomalous diffusion (class II) 


see fig. 4b 


see fig. lib 


Triangular 


Cl 7^ Cr — > No-diffusion (class IV) 


Cl = Cr — > one critical line — )• Anomalous diffusion (class II) 


see fig. 16b 
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